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ABSTRACT 

Context, Despite of the importance of magnetic fields for the full understanding of the properties of accreting Herbig Ae/Be stars, 
these fields have scarcely been studied over the rotation cycle until now. One reason for the paucity of such observations is the lack of 
knowledge of their rotation periods. The sharp-lined young Herbig Ae star HD 101412 with a strong surface magnetic field became 
I . in the last years one of the most studied targets among the Herbig Ae/Be stars. 

Aims. A few months ago we obtained multi-epoch polarimetric spectra of this star with FORS 2 to search for a rotation period and to 
constrain the geometry of the magnetic field. 

Methods. We measured longitudinal magnetic fields on 13 different epochs distributed over 62 days. These new measurements to- 
gether with our previous measurements of the magnetic field in this star were combined with available photometric observations to 
determine the rotation period. 

Results. The search of the rotation period resulted in P = 42.076 ± 0.017 d. According to near-infrared imaging studies the star is 
observed nearly edge-on. The star exhibits a single-wave variation of the longitudinal magnetic field during the stellar rotation cycle. 
These observations are usually considered as evidence for a dominant dipolar contribution to the magnetic field topology. 
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,— I \ 1. Introduction inclinations and vsinz, and from the knowledge of v eq values 

and stellar radii the rotation periods can be deduced. Our esti- 

. A number of Herbig Ae stars and classical TTauri stars are sur- mation of rotation periods of 1? Herbig Ae stars and four de _ 

L; ; rounded by active accretion disks and, probably, most of the brifj disk stars showed that for a major part of the studied sam . 
. ^ , excess emission seen at various wavelength regions can be at- ple the periods are of the order of a few days or fractions of 
^ ■ tributed to the interaction of the disk with a magnetically active days> and only for one Herbig Ae st ^ HD 101412> the expected 
»H i star (e.g. Muzerolle et al.|2004j. This interaction is generally re- riod was longer than 17 days (Hubrig et al 2QQ9) Among 
«? , ferred to as magnetosphenc accretion. Recent magnetosphenc the studied sampl6i HD 101412 also showed me largest lon . 
accretion models for these stars assume a dipolar magnetic field gitudinal magne tic field, (B z ) = -454+42 G, measured on low- 
geometry and accreting gas from a circumstellar disk falling bal- resolution p i ar i me tric spectra obtained with FORS 1 (FOcal 
hstically along the field lines onto the stellar surface. Reducer low-dispersion Spectrograph) mounted on the 8-m 
Despite of the importance of magnetic fields for the full un- Kueyen (UT2) telescope of the VLT. The subsequent spectro- 
derstanding of the properties of accreting Herbig Ae/Be stars, scopic study of twelve UVES and HARPS spectra of HD 101412 
these fields have scarcely been studied over the rotation cycle reve aled the presence of resolved magnetically split lines indi- 
until now. One reason for the paucity of such observations is the cating a variable magnetic field modulus changing from 2.5 to 
lack of knowledge of their rotation periods. The rotational pe- 3 5 kG ( Hu brig et al. l20T0l The presence of a rather strong mag- 
riod has not been known for any Herbig Ae star or debris disk netic field on the surface f HD 101412 makes it a prime candi- 
star. The v eq values can be estimated for stars with known disk date for studies of ±s impact of me mag netic field on the phys- 

ical processes occurring during stellar formation. In this work 

* Based on observations obtained at the European Southern we present for the first time a mean longitudinal magnetic field 
Observatory (ESO programme 085.C-0137(A)). 
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measurement series, obtained with the multi-mode instrument 
FORS 2 at the VLT. The data and available photometric observa- 
tions are used to obtain the rotation period and to put constraints 
on the magnetic field geometry. This is the first, necessary step 
for future more detailed magnetic studies of this remarkable star. 



2. Period determination 

Multi-epoch series of polarimetric spectra of the Herbig Ae star 
HD 101412 were obtained with FORS 20 on Antu (UT1) from 
2010 March 30 to 2010 June 1 in service mode. Using a slit 
width of 0V4 the achieved spectral resolving power of FORS 2 
obtained with the GRISM 600B was about 2000. A detailed de- 
scription of the assessment of the longitudinal magnetic field 
measurements using FORS 2 is presented in our previous papers 
(e.g., Hubrig et al. 2004a, 2004b, and references therein). The 
mean longitudinal magnetic field, (B z ), was derived using 
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where V is the Stokes parameter which measures the circular 
polarisation, / is the intensity in the unpolarised spectrum, g e g 
is the effective Lande factor, e is the electron charge, A is the 
wavelength, m e the electron mass, c the speed of light, dl/dA 
is the derivative of Stokes I, and (B z ) is the mean longitudinal 
magnetic field. 

An initial frequency analysis was performed on the longi- 
tudinal magnetic field measurements using a non-linear least- 
squares fit of the multiple harmonics utilizing the Levenberg- 
Marquardt method (Press et al. |1992| l with an optional possibility 
of pre-whitening of the trial harmonics. To detect the most prob- 
able period we calculated the frequency spectrum for the same 
harmonic with a number of trial frequencies by solving the lin- 
ear least-squares problem. At each trial frequency we performed 
a statistical test of the null hypothesis for the absence of peri- 
odicity (Seber [19771 1, i.e. testing that all harmonic amplitudes 
are at zero. The resulting amplitude spectrum clearly showed a 
dominant peak with an equivalent period of 42.098 d. 

To increase the accuracy of the rotation period determina- 
tion and to exclude the presence of other periods, we used pho- 
tometric data that include 920 V-band and 539 /-band obser- 
vations from ASAS (Pojmanski 2002) and 1426 /?-band obser- 
vations from the robotic telescope Pi of the Sky (Burd et al. 
2005, Malek et al. 2010). Further, 78 triads of measurements 
(UBV bands) over eight nights with internal accuracy of 9.4, 
6.4, and 6.6 mmag were obtained in South Africa Astrophysical 
Observatory (SAAO) in April 2010. Combining the magnetic 
and photometric data, and using the method of Mikulasek et al. 
( 120 101 ). we obtain 



<y&/> n 



: MJD52797.4 ± 0.8 + 42.076 ± 0.017£ 



(2) 



The corresponding periodogram and the light curve varia- 
tions in the U, B, V, R, and / bands are presented in Figs. Q] 
and |2] The full description of the photometric analysis is given 
in Appendices A and B. A very asymmetrical shape of the R- 
band light curve, with a deeper minimum at phases 0.15-0.30 
and large scatter at the phase 0.85 is probably related to a specific 
wavelength range covered by this band, which includes the Ha 
emission line. As we showed in our previous work, this line is 




Fig. 1. Periodogram of HD 101412 built from VI ASAS data, 
data from Pi of the Sky, SAAO UBV, and (B z ) data. The 
periodogram displays only one prominent frequency / = 
0.0237665 d -1 , corresponding to the stellar rotational period of 
42.076 days. 
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1 The spectropolarimetric capabilities of FORS 1 were moved to 
FORS 2 in 2009. 
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Fig. 2. UBVR, and / light curves of HD 101412. Solid lines 
present the simplest periodic functions representing the observed 
photometric variations (see Appendix A). Because of the large 
scatter of individual measurements (see Appendix B) the stan- 
dard technique of normal points (as averages of several tens of 
phase adjacent measurements) was used. The light curves are 
shifted along the y-axis for clarity. 



doubled-peaked and strongly variable (see e.g. Fig. 1 in Hubrig 
et al. 120 101 ). The radial velocity of the blue component changes 
from -33.6 to -67.5 km s~' whereas the red component, which 
actually appears triple in the phase 0.17, is shifted by 25.5 to 
119.6 km s _I . Our measurements of the Ha emission line flux 
over the rotation cycle on available HARPS and UVES spectra 
presented in Fig.[3]do not contradict this proposition. However, 
the secondary minimum in the emission flux appears at phase 
0.74. The phases of the minima in the measured values roughly 
coincide with the phases where the longitudinal field approaches 
zero value, i.e. close to the magnetic equator. 

The logbook of the FORS 2 spectropolarimetric observations 
is presented in Table Q] In the first column we show the MJD 
values for the middle of each spectropolarimetric observation. 
The phases of measurements of the magnetic field are listed in 
Column 2. In Columns 3 and 4 we present the longitudinal mag- 
netic field (S z )an using the whole spectrum and the longitudinal 
magnetic field (B z ) hyd using only the hydrogen lines. In the first 
two lines of the table we list the earlier measurements published 
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Fig. 3. Variations of the Ha emission flux over the rotation pe- 
riod (open circles), compared to the variation of the /?-band light 
curve (filled circles). The internal errors of the Halpha emission 
flux measurements are ordinarily smaller that the symbol size, 
but gross variations of the flux from cycle to cycle are expected. 

Table 1. Magnetic field measurements of HD 101412 with 
FORS 2. Phases are calculated according to the ephemeris of 
MJD = 52797.4 + 42.076 E. All quoted errors are l<x uncer- 
tainties. 



MJD 


Phase 
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[G] 


54609.190 


0.060 


-312+32 


-454±42 


54610.081 


0.081 


-235±28 


-317+35 


55285.141 


0.125 


-132±54 


-149+68 


55286.115 


0.148 


25±48 


209+74 


55287.107 


0.172 


49±36 


96+61 


55311.142 


0.743 


-202±44 


-235+56 


55312.068 


0.765 
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-241+60 


55320.056 


0.955 
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-784+57 


55322.149 


0.004 


-360±50 
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0.030 


-393+51 


-515+75 


55324.073 


0.050 


-359+63 


-453+64 


55327.159 


0.124 


-25±49 


-44+59 


55334.071 


0.288 
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338+66 


55335.152 


0.314 


351+75 


389+84 


55348.097 


0.621 


168±44 


322+60 
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Fig. 4. Phase diagram with the best sinusoidal fit for the longitu- 
dinal magnetic field measurements using all lines (filled squares) 
and hydrogen lines (open circles). 
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Fig. 5. Variations of equivalent width, radial velocity, line width, 
line asymmetry, mean magnetic field modulus, and mean longi- 
tudinal magnetic field as a function of the rotational phase. 



by Hubrig et al. ([7009). The corresponding phase diagram for 
all available mean longitudinal magnetic field measurements us- 
ing the whole spectrum and for those using only the hydrogen 
lines, including the best sinusoidal fit, is shown in Fig. [4] The 
variation has a mean of (B z ) = 9 + 18 G and an amplitude of 



= 465 + 27 G. 



In Fig. [5] we present the variations of equivalent width, ra- 
dial velocity, line width, and line asymmetry of a sample of iron 
lines, the mean magnetic field modulus, and the mean longitu- 
dinal magnetic field as a function of the rotational phase. Due 
to imperfect phase coverage the interpretation of the observed 
variations of the iron line profiles is not straightforward, but it 
seems that equivalent width and line asymmetry minima appear 
in the phase close to the surface field maximum of 3.5 kG and the 
positive extremum of the longitudinal field, while the presence 
of a secondary minimum is likely at the phase of the negative 
extremum of the longitudinal field. The line width is the largest 
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at the phase close to the field positive extremum. The radial ve- 
locity variations are most difficult to interpret: it is possible that 
they show double modulation, but it is not well-resolved in our 
measurements due to insufficient phase coverage. 

3. Discussion 

HD 101412 exhibits a single-wave variation of the longitudinal 
magnetic field during the stellar rotation cycle. These observa- 
tions are usually considered as evidence for a dominant dipolar 
contribution to the magnetic field topology. Our recent study of 
the abundances of HD 101412 using UVES an d HAR PS spectra 
resulted in v sin i — 3 + 1 km s _1 (Cowley et al. 1201 Oi l. The incli- 
nation angles of disks of Herbig Ae/Be stars (which are expected 
to be identical with the inclination angle of the stellar rotation 
axis) can be reliably derived only for resolved observations of 
disks. Fedele et al. ( 2008) used VLTI/MIDI observations to de- 
termine /=80±7°. Using this value we obtain v eq =3±l kms~'. 
Assuming that HD 101412 is an oblique dipole rotator, we fol- 
low the definition of Preston ( 1967): 



(B z r 



cos B cos i - sin B sin ; 



<B z ) max cos B cos i + sin B sin i 
so that the obliquity angle B is given by 
1 



B = arctan 



l + r 



cot i 



(3) 



(4) 



From the phase curve determined above, with 
< J g z ) max =474+32G and <B z ) min =-456+32G, we find 
r=-0.962+0.075, which for /=80+7° leads to a magnetic 
obliquity of £=84+13°. 

We note, however, that there are some gaps in the phase 
coverage for the mean longitudinal magnetic field observations, 
especially the extrema of the variation curve are not very well 
constrained. The true magnetic field geometry can probably be 
resolved once additional longitudinal magnetic field measure- 
ments as well as magnetic field modulus measurements over 
the full stellar rotation cycle become available in future obser- 
vations. 

For the first time rotation modulated longitudinal magnetic 
field measurements and photometric observations were used to 
determine the rotation period in a Herbig Ae star. The detec- 
tion of a large-scale organised predominately dipolar magnetic 
field on the surface of the young Herbig Ae star HD 101412 con- 
firms the scenario that the accretion phenomenon in young stel- 
lar objects is magnetospheric accretion. However, we are still at 
the very beginning of learning how magnetic fields of Herbig 
Ae stars can be incorporated in the modeling of stellar magne- 
tospheres which takes into account the full complexity of the 
circumstellar environment, including the observed outflows and 
jets in these objects. 

Longitudinal magnetic fields of the order of a few hundred 
Gauss have been detected in about a dozen Herbig Ae stars (e.g., 
Hubr ig et al . l2004cl [20061 [20071 [20091 Wade et al . 120071 Catala 
et al. 120071 1. For the majority of these stars rather small fields 
were measured, of the order of only 100 G or less. Due to the 
presence of the rather strong magnetic field in the atmosphere 
of HD 101412, and a very low vsin/ value of 3±lkms _1 , this 
star is one of the most suitable targets to study its environment, 
which includes the magnetosphere, the accretion disk, and the 
disk wind, altogether producing prominent emission features in 
the hydrogen line profiles. 
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4. Appendix A: The search for periodic variations 

The continuously updated archive of the All Sky Automated 
Survey (ASAS) is a useful source of photometric observations. 
Observations were obtained from two ASAS observing stations, 
one in LCO, Chile (since 1997) and the other on Haleakala, Maui 
(since 2006). Both are equipped with two wide-field instruments, 
observing simultaneously in V and / band. More details and 
the data archive are on http://www.astrouw.edu.pl/asas/ After 
removing observations of declared low quality (C, D) and ap- 
parent outliers, we obtained a set of 539 /-band and 920 V-band 
measurements well covering the periods 1997-1999 and 2000- 
2009, respectively. The robotic telescope Pi of the Sky has been 
designed to monitor a significant fraction of the sky with good 
time resolution. The final detector consists of two sets of 16 cam- 
eras, one camera covering a field of view of 20° x 20°. The set of 
the HD 101412 measurements covers the time interval of 2006- 
2009. More details can be found at http://g rb7uw.edu.pl/ CCD 
photometry was done without any filter, so that the results are 
similar to a broad-band R colour. The star was also observed over 
eight nights in 2010 April employing the 0.5 m reflector with the 
classical photometer in SAAO in UBV, using a fairly conven- 
tional single channel photometer with a Hamamatsu R943-02 
GaAs tube. We obtained 78 triads of measurements with an in- 
ner accuracy of 9.4, 6.4, and 6.6 mmag, respectively. 

The very good initial estimate of the period allows us to 
describe the observed periodic variations by a series of phe- 
nomenological models described by a minimum number of free 
parameters, including the period P and the origin of epoch count- 
ing Mo. The behaviour of the light curves in V and / is nearly the 
same. For simplicity we assumed a similar behaviour also for 
light curves in U and B bands, while the light curve R behaves 
differently. The periodic component in light variations can then 
be described by means of periodic functions F( ff) and Fr($r): 
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The functions F(§) and F r ($r) are the simplest normalised 
periodic function that represent the observed photometric vari- 
ations in detail. The phase of the brightness extreme is defined 
to be 0.0, and the effective amplitude is defined to be 1.0. The 
functions, being the sum of three terms, are described by two di- 
mensionless parameters B\ , Bo and B 3 , B\. § and § R are the phase 
function. Assuming linear ephemeris, we respectively obtain: 
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Fig. 6. Long-term changes in AS AS V (open circles) and / (full 
squares) corrected for variations. 



= (f - M () ) IP, & R = J?- A/*, (7) 

where A/r is the phase shift of the basic minimum of the func- 
tion Fr versus zero phase. 

Periodic changes of magnitudes in U, B, V, I: rrij(t) and 
changes in R: m R {t) are given by the relations: 

mj {t) = m-+AF{ff), m R {t)=rni + AF{ff)+A R F R {§ R ), (8) 

where A is the semiamplitude of light changes common to all 
bands and A R is the amplitude of an additional component of 
light variability being non-zero only in R. mj and Tmr are mean 
magnitudes in the individual bands. 

The periodic variations of the mean value of the longitudinal 
component of the magnetic field intensity (B z ) derived from all 
lines can be well approximated by the simple cosinusoid: 

(B z ) = {B Z )+A m cos [2n (tf - Af m )] , (9) 

where (B z ) is the mean magnetic field intensity, A m is the semi- 
amplitude of the variations, and Af m is the phase of the magnetic 
field minimum. 

All 1 8 model parameters were computed simultaneously by 
a weighted non-linear least-square method regression applied to 
the complete observational material representing in total 3134 
individual measurements. We found the refined value of the ro- 
tational period P = 42 ( ?076(1 7) and the origin of phase counting 
put at the UBVI light maximum M Q = 2452797.9(8). 

JD max = 2 452 797.9(8) + 42<?076(17)£. (10) 

Parameters describing the functions F and F R are: B\ = 
0.29(13), B 2 = -0.49(10), B 3 = 0.36(16), andyS 4 = 0.02(18), the 
function F R primary minimum phase Af R \ - 0.24(4), the phase 
of the secondary minimum A/r2 = 0.72(12), the semiamplitude 
of light changes A = -4.7(5) mmag, the semiamplitude A R = 
14.1 (2.2) mmag, V = 9.2660(14)mag, B = 9.4501 (13) mag, and 
U = 9.6106(18)mag. The phase of the minimum of the mean 
projected intensity of the magnetic field (B z ), Af m = -0.091(33), 
the mean value of it: (B z ) - 9 ± 18G, and the semiamplitude 
A m = 465 ± 27 G. 



5. Appendix B: Non-periodic variations 

Relatively precise ASAS V and / data show an apparent evolu- 
tion of the mean values to the extent of several hundredths of a 
magnitude over time scales of several years (see Fig.|6|. For this 
reason, we use detrended magnitudes, i.e. magnitudes where we 
have removed the long term trends, in our search of periodic 
stellar variability. 

An inspection of the light variation in R plotted over the rota- 
tional phase in 2006-2009 reveals dramatic seasonal variations 
in the shape of the light curve. Unfortunately, the low accuracy 
of the R photometric measurements together with the observed 
long-term instability and the large stochastic changes do not al- 
low us to analyse this phenomenon more deeply. Using principal 
component analysis tools, we tested the possibility of seasonal 
changes in the shape of the V light curve in the first decade of 
this century. We conclude that the V light curve as defined by 
ASAS and SAAO measurements was quite stable during this pe- 
riod. 

Herbig Ae/Be stars are known to vary in a complex way. 
The mechanism causing this variability is not well understood. 
These more or less stochastic variations amount up to several 
hundreds of a magnitude on time scales from several minutes to 
tens of years (see e.g. Ruczinski et al. 2010). The light changes 
of HD 101412 in V and / presented in Fig. [6] can be treated as 
stochastic variations on the scale of tens of years. Nevertheless, 
there is also evidence of such changes on much shorter time- 
scales. 

Our own 78 UBV photoelectric measurements in 2010 April 
in SAAO were obtained in a total of 12 sets over eight nights. 
The duration of one measurement set was not longer than 15 
minutes and we assumed that the brightness of the star did not 
change during the exposure. Then we can relatively well esti- 
mate the typical inner uncertainty of one individual measure- 
ment as 9.4, 6.4, and 6.6 mmag in U, B, and V, respectively. 
However, the means in particular sets in U, B, and V, consist- 
ing of six individual measurements in the average, show an ad- 
ditional scatter of 14, 10, and 1 1 mmag (see Fig.|7). We inferred 
that this scatter is due to stochastic changes on the scale of hours 
or tens of hours. It should be noted that the observed stochas- 
tic variations in particular colours are highly correlated, whereas 
their amplitude is the largest in the f/-band. Such a behaviour 
is certainly somehow related to the physics of the mechanisms 
causing such changes. 
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Fig. 7. Stochastic changes in our SAAO UBV observations. Note 
that the scatter in individual bands is highly correlated. 



Because of the circumstance that we do not have at our dis- 
posal measurements in V and / during one night, we estimate 
the inner accuracy of these measurements implicitly. The mean 
scatter of detrended ASAS data in V and / are 16 and 10 mmag. 
Assuming that the estimate of stochastic noise in V of 1 1 mmag 
is valid or smaller for all ASAS V measurements, we conclude, 
that the inner uncertainty of one ASAS V observation is smaller 
than 12 mmag. The ASAS estimates of uncertainty of / obser- 
vations should then be smaller by a factor of 0.6 compared to V 
observations. The inner accuracy should then be about 7 mmag, 
as well as the measure of the stochastic scatter. In all cases the 
stochastic scatter in / is not allowed to be larger than 10 mmag. It 
means that the tendency of the decrease of the stochastic scatter 
with the increase of the wavelength is supported also by ASAS 
data. A different amount of stochastic variation is expected in R 
band, where an additional variability mechanism is probably ac- 
tive. The mean scatter in the /?-band is about 56 mmag. The inner 
uncertainty of Pi measurements can be derived from magnitudes 
obtained in individual nights - it is estimated to be 49 mmag. 
The stochastic scatter then can be evaluated as 30 mmag. 



